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The halo nuclei 6 He and 8 He are described in a consistent way in a microscopic multiconfiguration model, the refined 
resonating group method. The ground state properties have been calculated, and momentum distributions of fragments and 
neutrons have been determined in a simple reaction scenario, taking into account final-state interactions. The correlation of 
neutrons and fragments are investigated. 
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I. INTRODUCTION 

Nuclear halos - a meanwhile well established phenom- 
ena of light neutron rich nuclei - are considered to be 
roughly understood; and at least the ground state struc- 
ture of some halo nuclei like 11 Be and especially 6 He are 
even thought to be almost completely enlightened. 

Also the measured momentum distributions of halo nu- 
clei fragments can be described very well within various 
theoretical models [1-7] . But in all these models some free 
parameters have been adjusted to reproduce the proper- 
ties of one nucleus. Contrary to this our aim is a consis- 
tent description of a large number of light nuclei, includ- 
ing halo nuclei, within one model whose parameters are 
fixed once and for all. 

First results of this project have been presented in ||] 
for a consistent description of the nuclei 4 He - 8 He. It 
turned out that the results for increasing mass number 
become more and more sensitive to the P-wave part of 
the employed effective nuclcon-nucleon interaction. This 
part of the potential leads to an overbound 8 He indicating 
a lack of repulsion and hence yielding a too small radius. 

In order to cure this problem we removed the odd par- 
tial wave part of the effective potential and added a re- 
alistic NN interaction with a more repulsive core. 

Here we show that with this modified interaction in 
addition to the lighter isotopes also the heavier isotopes 
can be described very well in a consistent way, without 
adjusting any parameter of the interaction to the indi- 
vidual systems. 

Besides the ground state properties we compare mo- 
mentum distributions of fragments and neutrons from 
6 He and 8 He breakup - calculated in sudden approxi- 
mation and taking into account final-state interactions - 
with experimental data and calculations in other models. 
We also discuss the correlation of the halo neutrons. 



II. SPATIAL STRUCTURE 

All calculations have been performed in a microscopic 
multiconfiguration cluster model - the refined resonating 
group method (RRGM). Details of the model are given 
in §. 

The 6 He wave function is a superposition of the 
(nonorthogonal) clusterings a(nn), (an)n and tt corre- 
sponding to the 4 He+2n, 5 He+n and 3 H+ 3 H thresholds 
of the 6 He system. On the cluster relative coordinates 
orbital angular momenta Z = 0, 1 have been allowed for 
the 6 He ground state. Between the nuclei inside the a 
and triton clusters only I = are considered. Four Gaus- 
sian functions were employed for the description of each 
cluster relative wave function and one Gaussian function 
for the a and the triton cluster each. All Gaussian width 
parameters have been determined by variation. 

TABLE I. Calculated binding energy _B t h with respect to 
the 4 He+2n threshold, matter r m and charge r c h r.m.s. point 
nucleon radii for 6 He compared with experimental data. 
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In Table | the calculated binding energy with respect 
to the first threshold and the matter and charge radii are 
compared with experimental data. The binding energy 
is about 0.1 MeV to small and hence the radii are a bit 
too large. 

Table |n] shows the remaining binding energy for the 
case that one omits each one of the configurations in or- 
der to determine the importance of each configuration. 
Without the (an)n clustering with P-waves on both clus- 
ter relative coordinates (the 5 He g.s. + n configuration) 
the binding is lost. Omitting the triton-triton clustering 
decreases the binding considerable, which is surprising 
because of the high threshold but in agreement with the 
calculations of Csoto || . Omitting any one of the other 
configurations effects the binding energy hardly at all. 
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TABLE II. Remaining binding energy for 6 He omitting 
each one of the configurations. 



omitted configuration 



E th [MeV] 



none 
(an)n 
(an)n 
tt 

a(nn) 
a(nn) 



h,2=0 
h,2=l 

1=0 

h,2=0 

h,2=l 



0.85 
0.81 
unbound 
0.22 
0.81 
0.67 



The P-state probability is 14.1 % and thus in agree- 
ment with a hyperspherical harmonics calculation jig ]. 

The 8 He wave function is a superposition of the 
(nonorthogonal) clusterings (a(nn))(nn) with I — on 
all coordinates, (((an)n)n)n with I = 1 on all coordi- 
nates, ((an)n)(nn) with I = 1 to the n clusters and / = 
to the nn cluster, ((a(nn))n)n with I = 1 to the n clus- 
ters and I — to the nn cluster, and ((tt)n)n with I = 
on all coordinates and in addition also I = between the 
t clusters and I = 1 to the n clusters. 

The calculated binding energy with respect to the 
4 He+4n threshold and the matter and charge radii are in 
good agreement with experimental data (cf. Tabic III). 



TABLE III. Calculated binding energy £th with respect to 
the 4 He+4n threshold, matter r m and charge r c h r.m.s. point 
nucleon radii for 8 He compared with experimental data. 
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/] r m [fm] 




r ch [fm] 


Experiment 3.112 




P 2.49±0.04 
2.45±0.07 


1< 


1 1.76±0.03 IIS 



RRGM 



2.99 



2.41 



1.71 



The most important configuration for the binding en- 
ergy is {{{an)n)n)n (cf. Table IV). The contribution of 
the {{tt)n)n channel is substantial. The P-state proba- 
bility is 49.8 %. We omitted the D-states totally since 
it turned out that they contribute only a few keV to the 
binding energy. 

Both wave functions are determined with the same NN 
interaction by a variational principle. For the even par- 
tial waves the effective potential of Stowe and Zahn p0| , 



TABLE IV. Remaining binding energy with respect to the 
4 He+4n threshold for 8 He omitting each one of the configu- 
rations. 



omitted configuration 



E th [MeV] 



none 

(a(nn))(nn) 
(((an)n)n)n 
((an)n)(nn) 
{(a(nn))n)n 
((tt)n)n 



2.99 
2.97 
unbound 
2.96 
2.96 
2.36 



for the odd partial waves the realistic interaction of Eike- 
meier and Hackenbroich fl~l| l has been employed. Both 
contain Coulomb, nuclear central, spin-orbit, and tensor 
components. The strength of the spin-orbit force was 
increased by a factor of 2.5. 



III. MOMENTUM DISTRIBUTIONS AND 
SPATIAL CORRELATIONS 

In the simple picture of the serber model (sudden ap- 
proximation) |13| in which any influence of the reaction 
mechanism of the fragmentation and the final-state inter- 
action is neglected, the measured momentum distribution 
is given by the Fourier transformation of the relative wave 
function between the fragments integrated over the not 
observed momentum coordinates: 



dN 
dPT 



V / dPydP z 



with 



— / dn*UR) exp(iRP) 



(i) 



(2) 



The totally antisymmetrized relative wave function be- 
tween the fragments \4>\) and \cf>2) is calculated as an 
overlap with the total wave function |<^ 3 ) where the in- 
tegration is performed over all coordinates except the 
relative coordinate r between the two fragments: 



<ZW(-R) = {<Pi<f> 2 \S(r - R)A[ 



(3) 



In this sense the relative wave function is projected out of 
the multiconfiguration total wave function. This is neces- 
sary since the different configurations are not orthogonal. 
A denotes the antisymmetrization operator. 

The relative wave function can be projected onto a 
certain orbital angular momentum L: 



S&(Jl) - Y LL {R) 



dR'Yl L {R')^{R'). 



(4) 



In the common reaction scenario one of the halo neu- 
trons is removed by the target and scattered to large 
angles and therefore not detected. This picture is based 
on the average neutron multiplicity measured to be close 
to one Jl6| . The remaining second neutron and the core 
may form a short living resonant state. The measured 
momentum distribution of the products of the decay of 
this resonance are thus no more only determined by their 
original, internal motion in the halo nucleus. Rather the 
final-state interaction has, at least in some cases, an im- 
portant influence. The modification of the internal mo- 
mentum distribution by the final-state interaction can be 
described ||(7 16 by using a Breit-Wigner formula, 



r 



2tt (E - P rcs )2 + ra/4 ' 



(5) 
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FIG. 1. Square of the norm of the projected relative wave 



function ||!f , r . 



as a function of the r.m.s. radius of Inn) 



with resonance energy E rcs and width r. 

For the calculation of the momentum distribution of 
the 4 He fragments from 6 He breakup one has \(f>\) — 
| 4 He), \(j, 2 ) = \nn) and |0 3 ) = | 6 Hc). The | 4 He) func- 
tion was determined by variation and has, in good agree- 
ment with experimental data, a point nucleon root mean 
square radius of 1.42 fm. The a cluster radii inside | 6 He) 
are very similar to this value but the momentum distri- 
bution is not very sensitive to it. 

To determine the function \nn) and thus the correla- 
tion of the two halo neutrons in 6 He we took the norm of 
the projected relative wave function ||^ re l|| as a measure 
of the content of a certain \nn) configuration in | 6 He) 
for a fixed | 4 He) fragment. Figure [l] shows H^ill 2 as a 
function of the root mean square radius of \nn) which 
is a good parametrization of the simple Gaussian func- 
tion which has been used for the r-space representation 
of \nn). The two pronounced peaks in Fig. l], at 1.0 fm 
for L rc i = and at 2.3 fm for L IC \ = 1, correspond to a 
spatialy strongly correlated dineutron configuration and 
a spatialy strongly anticorrelated cigar- like neutron pair. 
This observation is in agreement with the results of a 
hyperspherical harmonics calculation (HH) |13] . 

The \nn) configurations with maximal |^ r ~|"[| have been 
choosen for the calculation of !? re ] (R) . The resulting mo- 
mentum distribution for 4 He fragments in sudden ap- 
proximation is compared in Fig. |2| w ith the measured 
transverse momentum distribution |14] and a calculation 
in a HH model Jl2| where also the approximation of the 
serber model was used. Both theoretical curves are in 
good agreement with the experiment. This indicates 
that, besides the fulfilled preconditions for the serber 
model, high energy, light target and a loosely bound pro- 
jectile, the 4 He n final-state interaction hardly effects the 
4 He momentum distribution. Also the 4 He target interac- 
tion seems to play no important role for the detected 4 He 
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FIG. 2. 4 He transverse momentum distribution from a 6 He 
breakup calculated in the RRGM with final state interactions 
(thick solid line) and in sudden approximation (dashed line) 
in comparison with experimental data jw| and the HH model 
[ fl2| in sudden approximation (thin solid line). 

fragments. 

Taking final-state interactions into account, the He 
transverse momentum distribution is given by 



oc 



APydPz 



dp'dp" p, HC{nn) (p) 



Zbw(P') S(P' + -P" + P), 



(6) 



where P" is the momentum of the removed neutron and 
P the 4 He momentum, both in the 6 He cm. system (cf. 
Fig. ||). P' is the relative momentum of the detected 
neutron and the 4 He fragment in the 5 He cm. system, 
connected with E from Eq. (0) by 



E = 



pl2 

2§m 



with the nuclear mass m. We used E rcs = 0.89 MeV 
and r = 0.6 MeV in Eq. (ph for the 5 He ground state 
resonance [jlTj. Indeed Fig. g shows that the final-state 
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4 He 




detected n 

FIG. 3. Definition of the momenta appearing in Eq. (t 
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FIG. 4. Neutron transverse momentum distribution from a 
6 He breakup calculated in the RRGM with final state interac- 
tions (thick solid line) and in sudden approximation (dashed 
line) in comparison with experimental data [jl5| and the HH 
model |f| in sudden approximation (thin solid line). 

interaction does almost not effect the momentum distri- 
bution of the 4 He fragments. 

In order to calculate the neutron momentum distri- 
bution one has to employ |0i) = | 5 He) and |02) = \n). 
The wave function for the unbound system | 5 He) was 
constructed in the same way as for the dineutron \nn). 
The maximum of ||^rci|| was determined with respect to 
the relative wave function between the a cluster and the 
neutron in | 5 He) for a fixed a. In this sense the mean 
value of the a core halo neutron radius inside 6 He was 
determined to 2.0 fm. The corresponding momentum 
distribution (cf. Fig. ^) in sudden approximation is much 
broader than the measured transverse neutron momen- 
tum distribution [fl5| but in agreement with a three-body 
calculation in the serber model ||] . 

Some different processes might be responsible for this 
discrepancy. The 4 He n final-state interaction is consid- 
ered to be the most dominant one Q. But moreover, 
in contrast to the 4 He fragments the detected neutrons 
might have different origins: they could be emitted from 
the neutron halo by a direct breakup or with a target 
interaction, they could be emitted from the 4 He core or 
even from the target itself || . 

Taking final-state interactions into account, the neu- 
tron transverse momentum distribution is given by 



dA^i 
dP T 



j dPydP Z j dP'dP>5 Hen (P") 

I bw (P')6(P'-P-Ip"), 



(8) 



where P is the momentum of the detected neutron. The 
meaning of P' and P" is the same as in Eq. || (cf. Fig. |^). 
Contrary to the 4 He fragment, the neutron momentum 
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FIG. 5. 6 He transverse momentum distribution from 8 He 
breakup calculated in the RRGM with final-state interactions 
(solid line) and in sudden approximation (dashed line) in com- 
parison with experimental data juj. 

distribution is considerably influenced by the final-state 
interaction (cf. Fig. ||). The distribution @ is in good 
agreement with the measured distribution. The remain- 
ing discrepancy beyond 60 M^L might be due to neglect- 
ing the neutron target interaction flqj. 

The momentum distribution for He fragments from 
8 He breakup was calculated in a way analogical to the 
determination of the 4 He distribution from 6 He breakup: 
The determination of maxima of showed that the 

second halo neutron pair has a strong correlated compo- 
nent with a root mean square radius of 1.2 fm for spin 
s = and a less correlated component with a radius of 
1.8 fm for spin s — 1. The momentum distribution for 
6 He fragments calculated in the RRGM model in sud- 
den approximation seems to be a bit to narrow but is 
still compatible with the measured transverse momen- 
tum distribution |l(| (cf. Fig. |J). 

Taking final-state interactions into account, the 6 He 
transverse momentum distribution is given by 



d/V FSI 

dP r . 



J dP y dP z J dP'dP" p6 He („„)(P) 

I BW (P')6(P' + ^P" + P), 



(9) 



where the momenta P, P' and P" are defined analogi- 
cally to Eq. (|). 

For the parameters E rcs and r in Eq. (|^) we employed 
not only the experimental values for the 7 He ground state 
resonance (0.44 MeV and 0.16 MeV) but also the re- 
sults of an S-matrix analysis of the phase shifts calculated 
in the RRGM (0.442 MeV and 0.114 MeV) in @. The 
resulting momentum distributions differ only within the 
line width in Fig. || and |^. It is shown in Fig. | that 
the final-state interaction increases the width of the 6 He 
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IV. SUMMARY 
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FIG. 6. Neutron radial momentum distribution from He 
breakup calculated in the RRGM with final-state interactions 
(thick solid line) and in sudden approximation (dashed line) 
in comparison with experimental data |Q and COSMA cal- 
culation |16| in sudden approximation (thin solid line). 



It was shown that the ground state properties of halo 
nuclei and the momentum distributions of fragments and 
neutrons from breakup reactions can be well reproduced 
in a consistent way in the microscopic description of 
the refined resonating group method. The momentum 
distributions of 4 He and 6 He fragments from 6 He and 
8 He breakup, respectively, are in good agreement with 
the measured distributions and almost not influenced by 
final-state interactions. On the contrary, the neutron mo- 
mentum distributions can only be explained by taking 
final-state interactions into account. 
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momentum distribution in such a way as to yield good 
agreement with the measured one. 

The neutron momentum distribution was calculated by 
setting \4>i) = | 7 Hc) and \<f>2) = \n). The maximum of 
||$rei|| was determined with respect to the relative wave 
function between the 6 Hc function and the neutron in 
7 He) for a fixed 6 He. The mean value of the 6 Hen ra- 
dius was determined to 2.3 fm. Contrary to Eq. ([IJ) the 
momentum distribution for neutrons from 8 He breakup 
has been calculated by integration only over one momen- 
tum coordinate in order to compare with experimental 
data @. The RRGM and the COSMA calculations @ 
in sudden approximation are in good agreement (cf. Fig. 
^J) with each other but much broader than the measured 
radial momentum distribution. This discrepancy is pre- 
sumable due to the same reasons as for the 6 He breakup. 

Taking final-state interactions into account, the neu- 
tron radial momentum distribution is given by 



dA^i 
P t dP r 



x 



j dP z J dP'dP" pr Hcn (P") 



Jbw(P') 5{P' P -P"). 



(10) 



The momenta P, P' and P" are defined analogically to 
Eq. (|). The extracted width of 66 MsY i s a little larger 
than the measured one of 54 (cf. Fig. ^J). 

The origin of this small deviation from the experimen- 
tal data is not obvious. The width of the distribution 
given by Eq. jlQ), however, depends very strongly on the 
position of the^He resonance but not on the resonance 
width. For example lowering the parameter E Tes in Eq. 
(||) by about 100 keV yields perfect agreement with the 
data (see also |lq]). 
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